X-ray and ~f-ray Astrophysics of Galactic Sources 



1 



Low Energy X-ray Emission from Young Isolated 
Neutron Stars 

M. RUDERMAN 1 ' 2 ' 3 

1 Department of Physics, Columbia University 

2 Columbia Astrophysics Laboratory, Columbia University 
3 Institute of Astronomy, Cambridge 



ABSTRACT. A young neutron star with large spin-down power is expected to be closely 
surrounded by an e^ pair plasma maintained by the conversion of 7-rays associated with the 
star's polar-cap and/or outer-gap accelerators. Cyclotron-resonance scattering by the e~ and 
e + within several radii of such neutron stars prevents direct observations of thermal X-rays from 
the stellar surface. Estimates are presented for the parameters of the Planck-like X-radiation 
which ultimately diffuses out through this region. Comparisons with observations, especially of 
apparent blackbody emission areas as a function of neutron star age, support the proposition 
that we are learning about a neutron star's magnetosphere rather than about its surface from 
observations of young neutron star thermal X-rays. 



1. Introduction 

A young, isolated, magnetized neutron star (NS) is expected to have several potentially 
important sources of low energy (~keV) X-ray emission (cf. Figure 1): 

a) Thermal X-rays from the whole stellar surface as the NS cools down after a very 
high temperature birth in a supernova explosion. 

b) Thermal X-rays from its much smaller but hotter polar caps, additionally heated 
by backnow to the stellar surface of extreme relativistic electrons or positrons. These 
leptons come from e ± producing accelerators on the open B-field-line bundles which 
connect surface polar-caps to the spinning NS's distant "light-cylinder." These accel- 
erators may be just above the NS surface ("polar-cap accelerators") or in the NS's 
outer-magnetosphere ( "outer-gap accelerators" ) . 

c) Synchrotron radiation from newly created pairs as the e~ and e + radiate away 
their initial momenta perpendicular to local B. 

The angular and spectral resolution achieved by the Chandra and XMM satellites 
allows discrimination of such NS radiation from that of the supernova remnant in which 
the youngest (age t <, 10 4 yrs) NSs are usually embedded, and also thermal components 
(a and b) from the power low component (c); cf. Figure 2. In Section 2 we consider rather 
dramatic differences between observed thermal X-ray emission from young isolated hot 
neutron stars and the expected emission if the NS (near) magnetospheres were essentially 
empty. Section 3 discusses consequences for this emission if the NS magnetosphere within 
several NS radii of the stellar surface has an pair density (n±) large enough to give 
multiple cyclotron-resonance scattering of NS thermal X-rays before they ultimately 
escape or scatter back to the star. Section 4 estimates the steady state n± and its 
consequences for observed NS thermal spectra. Section 5 summarizes some conclusions 
and further problems. 

2. Observations and models 

Reported thermal X-ray luminosities (L^) of young isolated neutron stars are shown 
as a function of NS "age" (t) in Figure 3. (The "age" is generally that of the SNE 
envelope around the NS when that can be estimated. For some of the youngest NSs, 
J2323, J0852, J1617, and J0821 spin-period P with spin-down-rate P are not known so 
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Fig. 1. An outer-gap (OG) accelerator, extending along the open field line bundle on both 
sides of its intersection with the null-surface where Q ■ B = 0. GeV 7-rays directly from this 
accelerator may produce pairs close to the neutron star. 10 2 MeV 7-rays, curvature radiated 
by TeV e + (e~) flowing out of this accelerator down to the star, will be a very strong source 
for such e ± . On the other side of the star a polar cap (pc) accelerator may also be a strong 
source for 7 — > e + + e~ . (Because e pairs from OG accelerators can quench PC accelerators 
and vice versa, it is unlikely that a particular bundle of open field-lines would sustain both.) 



the canonical spin-down age (P/2P) cannot be used.) The computed thermal is that 
for the simplest, initially very hot, 1.4M neutron star model: a core consisting solely 
of neutrons, protons and extreme relativistic electrons, surrounded by a relatively thin 
more poorly conducting crust which keeps the NS surface temperature (T) very much 
less than that at the crust-core boundary (T c ). Although model T depends upon assumed 
NS surface detail, L\ from core-cooling should not as long as T c ^> T. Disagreement 
between model estimates and observations is not great, and usually compatible with 
uncertainties in the distance to the NS. (PSR 1055, for example, may well be half as 
far away as the value used in Figure 3, decreasing its L\ by a factor four. Some of 
the extremely young NSs may have apparent L\ which are strongly diminished when 
they are observed along particular lines-of-sight, as discussed in Section 5e.) Comparison 
between predictions and observations of Figure 3 do not discourage attribution of those 
Lx to the expected NS cooling. 

However, there are two great surprises in the observed thermal spectra of the NSs. 

1) None of the spectral features calculated for NS surface compositions other than 
H (and, presumably, He) are observed (with the remarkable exception of 1E1207 which 
may have a very different explanation to be discussed elsewhere). Indeed, a blackbody 
spectrum is generally an acceptable fit to observations and radiation parameters (surface 
temperature and radiating area) have usually been calculated for it. 

2) The effective blackbody radiating area (Abb) has the magnitude and age depen- 
dence indicated in Figure 4, hugely different from a simple age independent Abb ~ 
A N s = 4:ir(R NS ) 2 ~ 47r(10 km) 2 ~ 10 13 cm. Instead, observations give 

A NS U0 4 yrsy V ; 

for t <; 10 5 yrs; and Abb/A^s ~ 0(1) only for older NSs. If, despite the support of 
Figure 3 for the observed being mainly powered by core cooling, it is assumed 
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Fig. 2. X-ray count flux from the Vela pulsar (after removal of nebular emission) . (After Mori 
et al. 2002.) 




33 



\0 



10" 



Fig. 3. Thermal X-ray luminosity (L x ) of young isolated canonically magnetized neutron stars 
as a function of stellar age (t) Observed L x are from summaries of XMM and Chandra results 
(Becker & Aschenbach 2002; Pavlov et al. 2002) supplemented with data from Gotthelf et al. 
(2002) for PSR B1706, McGowan et al. (2003) for J0538, Mori et al. (2002) for Vela, and Wang 
& Halpern (1997) for Geminga. The solid line model curve for L x to which these data are 
compared is based upon Page (1998); cf. Tsuruta (1997). 
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Fig. 4. Blackbody fits for thermal radiation areas of isolated neutron stars (Abb), in units 
of a nominal neutron star surface area Ans = 47r(10km) 2 , as a function of neutron star age. 
(References for these data are the same as those for Figure 3.) The decreasing dash-dot diagonal 
is the model area of a hot polar cap of a central dipole. The upper dashed line is from the 
r = lCP 3 s e outflow model of Section 4; the lower solid line is from the r = s model in 
which e flow before e + + e~ -^7 + 7 annihilation is assumed to be ignorable. 



that the dominant sources of L 1 ^ are the hotter polar caps, the difference between 
expectations and observed black body areas becomes even greater: NSs with central 
dipolcs (or uniform interior magnetization) have polar caps of area (A pc ) with 



A pc 5 • 10" 



Ans P(sec) 



(2) 



This ratio is shown as the downward sloping dashed line of Figure 4, differing even more 
dramatically in both slope and magnitude from observations than does Abb/Ans ~ 1- 
If it is assumed that the NS surface, instead being a perfect absorber, has an H- 
atmosphere, L\ is not much changed. The associated affective radiating areas for near 
keV X-rays still grow almost linearly with t but now would have very significantly greater 
magnitudes. Typically the two different inferred NS surface temperatures, Tbb and Th, 
differ by about a factor two; Tbb = (1-5— 3)Th (Pavlov et al. 2002; Becker & Aschenbach 
2002) and have area ratios inversely proportional to the surface T 4 . Then for putative 
H-atmospheres when only blackbody parameters have so far been reported it seems 
plausible to estimate Ah ~ 2 4 Abb- the youngest NSs would still have Ah /Ans *C 1. 
The dramatic rise and subsequent leveling off in radiating areas would still remain and 
require explanation. 
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3. Thermal X-rays observed from a hot strongly magnetized neutron star 
embedded in a dense e ± plasma 

A surrounding e plasma with number density n± = n+ = ra_ < 10 17 cm~ 3 , extending 
out to several stellar radii around a NS, would not give important Thomson cross section 
scattering of thermal X-rays emitted from the NS surface. However, because the thermal 
X-ray's energy (Ex = huj ~ keV) is less than the e - and e+ cyclotron-resonance energy 
at the stellar surface (Tlujb = eBh/mc ~ lOB^kcV), an emitted thermal X-ray will 
pass through a cyclotron-resonance where Ex ~ ficosir) within a few stellar radii of 
the emitting surface (Wang et al. 1998, Zhu & Ruderman 1997,Dcrmer & Sterner 1991, 
Mitrofanov & Pavlov 1982). The resonant polarization-averaged X-ray scattering cross 
section there is well approximated as 

47r 2 e 2 

a ± —8[u-u B {r)\. (3) 

mc 

We assume below that equation (3) also adequately approximateds the scattering cross 
section for any polarized x-ray if the averaging is over its e~ and e + scattering cross 
sections. (This fails when the photon's momentum is almost perpendicular to local B 
together with its electric polarization being almost parallel to B. However, when the 
optical depth for such resonance scattering is very large, as in the main applications 
discussed below, or the near-star B-field is very irregular, the limit to accuracy from 
the use of Equation 3 does not substantively alter the conclusions which will be based 
upon it.) 

With the a of Equation 3, the optical depth for cyclotron-resonance scattering (OD) 

of an energy Ex X-ray is about (\x/^s) ~ 10 5 times that for non-resonant Thomson 
scattering: 

OD = / <j±n±dr (4) 

Jr ns 6mcE x 

as long as re, the distance from the NS center where Hujb = Ex, satisfies r B > Rns- 
(The denominator factor 3 is based on i3-field r-dependence of a central dipole.) 

The cyclotron-resonance backscatter radius t\b ~ 3Rns for an Ex ~ IkeV X-ray 
emitted from a NS with surface dipole field 3 • 10 12 G. This backscattering region around 
r ~ 3Rns forms a reflecting, slightly leaky Hohlraum container holding in radiation 
from the NS. Its X-ray transmission coefficient 

'"-GhW' (5) 

and the container X-ray reflection coefficient f = 1 — t ~ OD(OD + 2) _1 . When t <C 
R%s/ r B ~ 10 _1 , an X-ray emitted from the NS surface will almost always be reflected by 
the Hohlraum container and then reabsorbed by the NS before escaping the Hohlraum. In 
this regime the .ffo/i^raMm-contained X-rays acquire a Planck-blackbody energy spectrum 
independent of the NS surface's emission and reflection properties. It is sensitive only 
to the (now higher) surface temperature of the NS. 

In a steady state the total outflow luminosity is unchanged by the plasma in which 
the NS is embedded. Equivalently, the container has a surface emissivity e « t = 1 — f 
and 

L x = A BB o SB T% B = ^r 2 B eas B T BB (6) 

(with <jsb the Stefan Boltzmann constant). 
Then from Equations (5) and (6) 

A BB e47rr| ~ (7) 
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as long as the computed Rns of Equation (7) gives an Abb A:irR 2 NS . In the opposite 
limit, as OD drops below (rs / Rns) 2 /2, Abb — > 47r_R^ s . If Abb varies greatly with Ex 
observed emission out from the Hohlraum 's container would have a spectral distribution 
significantly different from the interior Planck one. 

4. Inferred radiation areas (Abb) and neutron star spin-down age 

In a steady state the production rate of e ± pairs near a spinning-down NS (N±) needed to 
maintain a constant local pair density n± in the relevant cyclotron-resonance scattering 
region is balanced by the sum of two local pair removal rates: 
a) Local annihilation into 7-rays (n±) with 



(assuming nonrelativistic e + — e relative velocities along B and anti-parallel magnetic 
moments along the strong local B have been achieved). 

b) Pair outflow (n£) from strong push of incident thermal X-rays on the continually 
replenished cyclotron-resonant container of the Hohlraum give a finite residence time 
(r) there: 



One estimate for the r follows from the assumption that very quickly after the pair 
creation the e^ 1 speed along local B becomes such that the push on e~ + e + is perpen- 
dicular to local B (about 10 10 cm s _1 ). Then to move along that B through the container 
thickness (~ tb ~ 3Rns ~ 3 • 10 6 cm) takes 



[Pair inflow from NS gravitational pull could be important if that pull were greater 
than the outward push by contained Hohlraum X-rays. In the models presented here the 
pressure for that outward push is about [OD /2)(L X / Airr B c) . The oppositely directed 
pull from the weight of the n± around t\b which constitute the container (Equation 4) 
is only 10~ 3 (10 33 erg s _1 / L x ) as large. Therefore we do not consider it further here. 
This interpretation does not support the gravity dependent e ± -layer structure and X- 
ray spectra proposed for the Crab pulsar by Zhu and Ruderman (1997); cf. also Wang 
et'al. (1998).] 

At the other extreme, if this local B near the NS is still so "complicated" that e ± is 
kept from flowing out before annihilation, 



The near-star pair input rate may be even more difficult to estimate than the an- 
nihilation rate. Pairs on closed field lines very near a young NS (r <; 3Rns) can be 
created by a) MeV - GcV 7-rays which enter into that region from outer-gap accelera- 
tors; b) MeV - 10 2 MeV curvature 7-rays from TeV e+(e~) flowing down to a NS polar 
cap on B-ficld lines connecting the polar cap to the outer-gap accelerators; and c) 1 - 
10 2 MeV 7-rays generated by polar cap accelerators on open-field lines, and then bent 
over to cross near-surface closed-field lines by the strong NS gravitational pull on them 
(a source which becomes large and important when the usual central dipole model for 
B is replaced by a more realistic one without axial symmetry). 

A further complication in calculating N± near the NS is a possible second generation 
of production by GeV 7-rays from outer-gap accelerators. The e 1 *" pairs they create, if 




(8) 



t - l(T 3 s. 



(10a) 



T ~ OO. 



(106) 
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Fig. 5. X-ray luminosity (Lx) for 2 keV < Ex < keV of canonically magnetized isolated 
neutron stars as a function of their ages (Possenti et al. 2002). The solid curve is Equation 12: 
L c = 10 33 (10 4 yrs/t) 2 erg s _1 . 



these 7-rays pass are within IORns of the NS, can themselves be the source of energetic 
(E~f > lMeV) synchrotron 7-rays. These, in turn, would create more pairs if they pass 
into the much stronger B much closer to the star (Cheng & Zhang 1999; Chang et 
al. 1998; Yakovlev et al. 2002). A rough estimate for the Vela pulsar suggests that 
N± ~ 10 37 - 10 38 s _1 for this pulsar. 

A more empirical, and probably more reliable, estimate of N± near the NS comes 
from identifying the Vela power law component of Figure 1 with the x-ray synchrotron 
radiation from the immediate quenching of e~ and e + perpendicular momenta of all 
newly created e 1 * 1 pairs in the strong local B (Wang et al. 1998). The power law syn- 
chrotron luminosity from the local field where the e are created by ~ 10 2 MeV 7-rays 
has a total power law luminosity (L X L ) for 2 keV < Ex < 10 keV 

/8kcV\ 

L? L (2 - lOkcV) ~ N±mc 2 ~ 10- 5 JV±erg s" 1 , (11) 

\ nu>B J 

where the local B field is assumed to be around 10 10 G so that %ujb ~ 10 _1 keV. If the 
L X L is indeed from such synchrotron radiation, its X-ray photon flux spectrum between 
2 and 10 keV should be ~ E x r with r ~ 1.5 as long as all relevant Tilob <, 2 keV and 
maximum synchrotron radiation energy Ex > 10 keV. [Although reported T ranges do 
not always include this value it is compatible with many examples: 1.25 - 1.45 for Vela, 
1.4 - 1.5 for PSR 0656, 1.6 for the Crab (Becker & Aschenbach 2002) and 1.66 for PSR 
1055 (Pavlov et al. (2002).] 

Observed (2 - 10 keV) Lx as a function of NS ages are shown in Figure 5. (In this 
range for Ex, Lx is expected to be dominated by its power law component.) 

We approximate a fit to these data by the line shown there: 

pl _ T __ _ ln 33 /^10 4 yrs 



L x ^ Lx = 10 3i ^^f^j ergs- 1 . (12) 

When combined it with the defining relations B\ = i?~ 6 c 3 0~ 4 x (spin-down power 
and NS spin-down age t = -0/20, Equation 12 becomes L x ~ ^BjR 6 NS I~ 2 x | 
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spin-down power|. For canonical pulsars with B^ ~ 10 12 G it is then roughly equivalent 
to the empirical Becker- Triimper relation (Becker & Triimper 1997; Seward & Wang 
1998; Possenti et al. 2002.) L x ~ 10~ 3 |spin-down power|. 

To calculate the needed n± by combining the total pair creation rate N± from Equa- 
tions 11 and 12 with the local disappearance rates of Equation (8, 9 and 10), it is 
necessary to estimate the volume (V±) in which the pair creation rate of N± is achieved. 
This depends upon the characteristic B in which the 7 + B — > e ± + B occurs, the same 
B as that of uib in Equation 11. For that Uujb ~ 10 _1 keV, corresponding to produc- 
tion beginning at B ~ 10 10 G by characteristic 10 2 MeV curvature 7-rays from the TeV 
e~(e + ) flow outward from a polar-cap accelerator and inward from an outer-gap one, 

V±(oc lu b 3 ) ~ 10 21 cm 3 . (13) 

Combining Equations 7 and 4 with re = 3 • 10 6 cm and Ex = lkcV, together with 
Equations 8, 9, 11, 12, and 13, gives n±, OD, and the two Abb/ Am s model curves shown 
in Figure 4. To obtain the solid curve we pretend that equation 10b is the appropriate 
entry for Equation 9, i.e., a very complex closed B-field structure to r ~ Tb- The opposite 
t = 10 _3 s approximation of Equation 10a gives the upper dotted line. The truth may 
well lie somewhere between the curves shown. 



5. Some conclusions, speculations, and problems 

The model results for estimated Abb magnitudes and their age dependence are qual- 
itatively similar to the observed ones of Figure 4 over several orders of magnitude. 
This gives considerable support to our central proposition. Strongly magnetized NSSs 
with large spin-down powers continually produce enough pairs on the closed B-field 
lines in their strongly magnetized near-magnctosphcres that X-ray cyclotron-resonant 
scattering there preempts direct observation of the NS surface. 

Although comparisons between these over-simplified models and NS thermal X-ray 
observations may not be quantitatively reliable they encourage further model-based 
interpretations of existing thermal X-ray observations. 

If the t = 10 _3 s model (uppermost curve) in Figure 4 leads to a reasonable approx- 
imation for cyclotron-resonance opacity around a young neutron star, that opacity 
alone would not explain why the NSs older than 10 4 yrs wold have L x fits to near 
blackbody spectra and Abb I Am s so much less than unity. These two features suggest 
that for those NSs (but not for younger ones) it may be more appropriate to fit observed 
thermal spectra with unobstructed emission from H or He atmospheres. (This seems, 
however, often to give optimal fits with Ah/Ams enough larger than unity to introduce 
new problems.) 

If the t — ► og model is a better description, surrounding e ± plasma blackbody-like 
emission from a slightly leaky e ± Hohlraum container would remain appropriate for very 
considerably older NSs (lower curves of Figure 4). However, at least one of these, J0538, 
would be expected to have such small surrounding opacity that the observations still 
suggest its interpretation as a NS emitting unobstructed thermal X-rays from an H (or 
He)-atmosphere with 

a h ( Tbb \ A Abb ^ Abb 1 ^ 
Am s \ Th J A N s A N s 

A more accurate description of the large n± around very young magnetized spinning- 
down NSs will be needed for many observational details. 

a) Variation of n± and thus OD, with (angular) location. This would give spin-phase 
modulation to the observed thermal flux. 
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b) The dependence of r 2 B /OD on X-ray frequency, causing some departure from a 
Planck spectrum in emission. (If the N± which sustains n± is inside rs and n± outflow 
through relevant tb conserves rirr 2 n±. Equation 4 does give mr B /OD independent of 
X-ray energy when B oc r~ 3 .) Another neglected departure of thermal X-ray emission 
from exact Planck spectra comes from significant e + /e~ flow velocities along B: e + /e~ 
recoil is very small in each cyclotron-resonance X-ray scatter (Ex/mc 2 <C 1) but not 
X-ray energy boost or degradation. 

c) Inclusion of backscattcr of hotter polar-cap X-rays which suppresses discrimination 
between the separate contributions to Lx of polar cap emission and general surface 
cooling. 

d) Correlations among L X L , Abb, and dipole moment for individual NSs. For ex- 
ample, special consideration must be given to relating N± to L X L from pulsars with 
exceptionally large 10 6 G) near their light cylinders. When such pulsars have strong 
outer-gap accelerators in that region, GeV (curvature) 7-rays + kcV X-rays produce 
pairs there with very strong initial synchrotron radiation in the 2-10 keV range. For 
these e^, local Tilob ~ 10 _2 eV so Equation 11 becomes N± ~ ltfL x L (erg s" 1 ) instead 
of the much larger inferred N± from L X L if the pairs are produced very near the NS. 
This distinction should be important for the Crab pulsar with its exceptionally large 
light-cylinder B (~ 10 6 G). The Crab pulsar's observed L P X L (2 - 10 keV), ~ 10 36 erg s" 1 , 
should and does lie well above the L cx t~ 2 line in Figure 5. Its L X L should not be used 
without adjustment to predict N± in the Crab pulsar magnetosphere and especially 
not in the inner-magnetosphere near the star. (A total crab pulsar N± inferred from 
Equation 11, without taking into account its large expected outer-magnetosphere con- 
tribution to L X L , would exceed observational upper limits for the Crab pulsar's half 
MeV annihilation 7-rays from e + + e~ -^7 + 7, (0.5 - 1) xl0 40 s" 1 (Ulmer et al. 2001; 
Massaro et al. 1991; Zhu & Ruderman 1997). 

e) Differences between total L x and that inferred from observations because the spin 
cycle averaged OD may vary greatly for observers along different lines-of-sight to the 
NS center. In particular, there is the possible escape of Hohlraum X-rays through an 
open field- line bundle "hole" (area ~ irrgilc^ 1 ) in the Hohlraum container. (Because 
of the large very relativistic particle flow in the bundle sustained by the accelerator 
somewhere on it, OD along it can be greatly suppressed at r ~ rs for thermal x-ray 
frequencies. Then, for pulsars with relatively large "hole" areas and large n± elsewhere, 
Hohlraum X-rays may mainly pass out through the container's polar-cap-like "holes" 
rather than diffuse out through the container's reflecting "wall" . When this obtains, the 
concentration of L x outflow through such holes could increase its apparent value by 
about a factor 2, or decrease it by a very much greater factor, depending on polar cap 
locations and that of the line-of-sight to the observer. (The rapidly spinning NS in 3C58 
(J0205 with P = 66 ms) with very large spin-down power (2 x 10 37 erg s" 1 ) would be 
expected to have a very large N± and thus n±. X-ray outflow through container holes 
should be greater than the diffusive one. To an observer for whom the container holes 
are largely hidden during the full NS spin cycle most of that NS's L x could be directed 
away from the observer. If so, the inferred L x magnitude would be greatly reduced from 
its true value. This may well account for not yet observing any L x from it (cf. Slane et 
al. 2002). 

Effects of N± may be complicated and difficult to evaluate quantitatively but their 
consequences may be crucial for interpreting observations of low energy X-ray fluxes 
from young isolated neutron stars. 
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